Quantum Coherence and Entanglement in the Avian Compass 



James A. Pauls, 1 Yiteng Zhang, 2 Gennady P. Berman, 3 and Sabre Kais 4, 5 E1 

1 Goshen College, Goshen, IN 46526 USA 
2 Department of Physics, Purdue University, West Lafayette, IN, 47907 USA 
3 Theoretical Division, T-4, MS B-213, Los Alamos National Laboratory, Los Alamos, NM 87545 USA 

4 Department of Chemistry, Department of Physics and Birck Nanotechnology Center, 

Purdue University, West Lafayette, IN 47907 USA 

5 Qatar Environment and Energy Research Institute, Qatar Foundation, Doha, Qatar 

The radical pair mechanism is one of two distinct mechanisms used to explain the navigation of 
birds in geomagnetic fields. However, little research has been done to explore the role of quantum 
entanglement in this mechanism. In this paper, we study the lifetime of radical pair entanglement 
corresponding to the magnitude and direction of magnetic fields to show that the entanglement 
lasts long enough in birds to be used for navigation. We also demonstrate that, due to a lack 
of orientational sensitivity of the entanglement in the geomagnetic field, the birds are not able 
to orient themselves by the mechanism based directly on radical-pair entanglement. To explore 
the entanglement mechanism further, we propose a model in which the hyperfine interactions are 
replaced by local magnetic fields of similar strength. The entanglement of the radical pair in this 
model lasts longer and displays an angular sensitivity in weak magnetic fields, both of these factors 
are not present in the previous models. 



Introduction. — The ability of many animal 
species, such as birds, insects, and mammals to 
sense the geomagnetic field for the purpose of 
orientation and navigation has led to huge in- 
terest in the field of biophysics [l], Q • There are 
currently two leading hypotheses to explain this 
remarkable ability: the magnetite-based mech- 
anism, and the radical pair mechanism 0-01 • 
The latter mechanism has been supported by 
results in the field of spin chemistry [H, [j| and 
by biological experiments [Icj. Recently sev- 
eral authors have raised the intriguing possibil- 
ity that living systems may use nontrivial quan- 
tum effects to optimize their orientation behav- 
ior IIEIEI. 

It has been suggested that entanglement, 
rather than mere quantum coherence is the con- 
tributing factor which allows the avian compass 
to achieve its high level of sensitivity Q . If this 
is so, does the duration of the entanglement last 
long enough to impact biological processes, and 
is the entanglement sensitive enough to the in- 
clination of the radical pair with respect to the 
Earth's magnetic field? To answer these ques- 
tions, we examine the lifetime of radical pair 
entanglement corresponding to different mag- 



netic field strengths, and compare the results 
with the candidate chemical reaction We 
also study the angular dependence of the radical 
pair entanglement within the geomagnetic field. 
Based on our results, we propose a new model 
to explore the underlying details. We find that 
the entanglement present in our proposed model 
displays both directional sensitivity as well as a 
sufficiently long duration of entanglement. 

Model. — The basic scheme of the Radical 
Pair Mechanism (RPM) involves three steps 
13 . 3]. The first step is light absorption, which 



is then followed by the formation of the radical 
pair and its interconversion between the singlet 
and triplet for electron spins. The final step 
is the decay of the singlet and triplet states 
to chemical products which produce a chemical 
signal detectable through an, as yet unknown, 
biological pathway. Typically, the radical-pair 
reaction involves two kinds of molecules which 
play the roles of electron donor and accep- 
tor. After absorbing light, the electronic state 
of the donor molecule is excited, e.g. one 
of the electrons in the donor molecule is ex- 
cited from the highest occupied molecular or- 
bital (HOMO) to the lowest unoccupied molec- 



ular orbital (LUMO) [14|. In the following step, 
if the two molecules are close enough to each 
other, the donor electron will be transferred to 
the acceptor molecule. Thus, both molecules 
will contribute an unpaired electron to form a 
pair of radicals. These two unpaired electrons 
in the donor and acceptor molecules are initially 
bound in a singlet state before they are spatially 
separated. Under the effect of the external mag- 
netic field, i.e. the geomagnetic field and local 
nuclear spins, the electron pair tranfers between 
the singlet and triplet states. Finally, the sin- 
glet and triplet states will produce different re- 
action products [H, 14 1 . The whole process is 
shown in Fig. [TJ 
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FIG. 1: Scheme of RPM. After a light-induced elec- 
tron excitation, the donor transfers an electron to 
the acceptor, forming a radical pair with the accep- 
tor molecule. The interconversion between singlet 
and triplet is affected by an external magnetic field. 
Finally the singlet and triplet decay into different 
products. 



To investigate the role of entanglement in this 
chemical compass, we focus on the second step 
of the RPM scheme, since it is while the elec- 
trons are separated that it is believed the ge- 
omagnetic field may influence intersystem con- 
version. 

Based on earlier literature jTij], we include 
only the Zeeman interaction and the hypcrfinc 
interaction in the Hamiltonian of the system 
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In Eq. (JTJ), the first term accounts for the Zee- 
man interaction, and the second term for the 
hyperfine interaction. (We assume that each 
electron is coupled to a single nucleus.) /; is 
the nuclear spin operator; Si is the electron spin 
operator, i.e., S = a/2 with a being the Pauli 
matrices; g is the p-factor of the electron, which 
is chosen to be g = 2^ ps is the Bohr magneton 
of the electron; and Ai is the hyperfine coupling 
tensor, a 3x3 matrix. 

As suggested in Ref. [Tj3|, we model the 
radical-pair dynamics with a Liouvillc equation, 



p(t)=--[H,p(t)] 
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In Eq. <j2j) , H is the Hamiltonian of the system; 
Q s is the singlet projection operator, i.e. Q s = 
\S)(S\, and Q T = \T + )(T+\ + \T )(T \ + \T__)(T-\ 
is the triplet projection operator, where \S) 
stands for the singlet state and (\T + ), \Tq), |T_)) 
stand for the triplet states [l^l; p(t) is the den- 
sity matrix for the system; k$ and kj- are the de- 
cay rates for the singlet state and triplet states, 
respectively. 

Calculations and Results. — For our calcu- 
lations we assume that the initial state of the 
radical pair is a perfect singlet state, |S) = 
TfGti) - lit))- Therefore, the initial condition 

for the density matrix is: p(0) = jIn <8> Q S , 
where the electron spins are in the singlet states, 
and nuclear spins are in a completely mixed 
state, which is a 4x4 identity matrix. As- 
suming that the recombination rate is inde- 
pendent of spin, the decay rates for the sin- 
glet and triplet should be the same [l3|, ks = 
kx = k, i.e., k is the recombination rate for 
both the singlet and triplet states. The ex- 
ternal weak magnetic field, B, representing the 
Earth's magnetic field in Eq. (jTJ) , depends on 
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the angles, 6 and ip, with respect to the ref- 
erence frame of the immobilized radical pair, 
i.e., B = Bo (sin 9 cos <p, sin 6 sin (p, cos 6>), where 
Bq = 0.5G is the magnitude of the local ge- 
omagnetic field. Without losing the essential 
physics, (p can be assumed to be 0. 

Since the radical pair must be very sensitive 
to different alignments of the magnetic field, it is 
necessary to assume that the hyperfine coupling 
tensors in Eq. ([1]) are anisotropic. However, for 
the sake of simplicity, we employ the hyperfine 
coupling as anisotropic for one radical, and the 
other as isotropic [13|, i.e., 
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Q T = \T)(T\, and |T) = \T+) + \T ) + |T_). 
The effect of the radical pair decay rates on the 
triplet yield has a twofold function [l3|. For a 
very high decay rate, i.e., larger than 10/zs -1 , 
the rapid decay of the radical pair prevents ef- 
ficient singlet-triplet mixing, as can be seen by 
the increase of the triplet yield in the weak mag- 
netic field. This means that the weak mag- 
netic field has very little effect on the triplet 
yields with fast decay rates. However, for very 
slow decay rates, i.e., smaller than 0.1/xs -1 , the 
triplet yield increases up to its maximum almost 
immediately when the magnetic field increases 
from zero, but remains essentially static as the 
magnetic field continues to increase. A decay 
rate of the order of l/^s -1 , seems to be opti- 
mum for the detection of a weak magnetic field. 
For all further calculations with this model we 
have assumed this value for our decay rate, i.e. 
k = 1/xs -1 . 
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FIG. 2: The triplet yields for three different decay 
rates as a function of the external magnetic field 
magnitude. The black line, for which k — 1/is -1 
seems to be a reasonable curve. 



In order to determine what values of the de- 
cay rates are reasonable for biological systems, 
we calculate the influence of different decay 
rates on the triplet yield, <I>t, as the exter- 
nal magnetic field varies, by using the parame- 
ters defined above. We define the triplet yield 
as [13 0: $t = k Tr[Q T p(t)]dt, where 
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FIG. 3: Entanglements for different magnitudes of 
the magnetic field for the angle 68° between the z 
axis of the radical pair and the magnetic field. 



Having fixed the decay rate to be 1/xs -1 , we 
study the radical pair entanglement as a func- 
tion of the magnitude of the geomagnetic field. 
The z axis of the radical pair is aligned at an an- 
gle of 68° with the magnetic field vector, which 
is the angle at which an earth-strength magnetic 



3 



field produces the largest triplet yield [13( . (See 
Fig. [4] in |13j.) In this paper, we use nega- 
tivity as the metric of entanglement, N(p) = 



II 1—1 



where ||p T ' 4 ||i is the trace norm of 
the partial transpose of the system's density 
matrix [l5j |. The results are shown in Fig. [3) 
We can see that when the magnetic fields are 
weaker than the Earth's magnetic field, or as 
strong as 1G, the entanglement curves are al- 
most identical. There does not appear to be 
any unique behavior that distinguishes a field 
in the neighborhood of 0.5 Gauss. However, 
under the Earth's magnetic field, the entangle- 
ment will be robust periodically during the first 
0.5/^s, which is longer than the suggested dura- 
tion of radical pair separation A stronger 
magnetic field (e.g. 5G) will disturb this peri- 
odicity. Previous research on the magnetic-field 
sensitivity of the chemical compass has demon- 
strated that the entanglement is helpful only if 
nature allows birds to optimize this behavior [H[ . 
On these grounds one can say that the entan- 
glement lasts long enough to play crucial role in 
the orientation of birds. 
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FIG. 4: Angular dependence of the triplet yields. 
The triplet yields are symmetric about 90°. 



So, we can say that the entanglement could 
play a role in the orientation and navigation of 
birds. We now recheck one of the properties of 



the avian compass, that it should depend on in- 
clination but not polarity. In Fig. [4] we see 
that the triplet yields are symmetric about 90°. 
Consequently, the radical pair mechanism can- 
not distinguish between magnetic fields that are 
oppositely directed but have the same magni- 
tude [H . 

The surprising result, in Fig. [5l is that the 
dynamics of entanglement does not change with 
angle, i.e., entanglement is not sensitive to the 
angle between the z-axis of the radical pair and 
the Earth's magnetic field. Therefore, it is rea- 
sonable to conclude that the entanglement of 
the radical pair cannot provide the same infor- 
mation as the triplet yields. In other words, 
using this model, the entanglement of the sys- 
tem does not directly affect the birds' ability 
to orient themselves. However, there might be 
indirect mechanisms which allow the birds to 
utilize entanglement. 
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FIG. 5: Entanglements for different angles. All 
curves are practically identical. In the geomagnetic 
field, entanglement does not change with orienta- 
tion. 



The above result (Fig. [5]) shows that the 
dynamics of entanglement are nearly static for 
different angles under the symmetric hypcrfinc 
tensors. This raises the question, what will 
happen when we use an asymmetric hyperfine 
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tensor? We examined several such cases. The 
asymmetric hyperfine tensors we applied are, 
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FIG. 6: Entanglements for different angles under 
the hyperfine coupling tensors A\. 

From Fig. [6] and Fig. [3 we can easily con- 
clude that the hyperfine coupling tensor pair of 
A? gives an intriguing result, the dynamics of 
the entanglement is clearly dependent on the 
system's orientation. This result inspired us to 
develope a new model in which only the exter- 
nal magnetic fields are considered, as in the case 
of A\ there are only two non-zero terms totally. 

New Model. — Previously, we had assumed 
that one electron of the radical pair experi- 
ences an anisotropic hyperfine coupling, while 
the other experiences an isotropic one. How- 
ever, this model cannot produce an angular- 
sensitive entanglement. On the other hand, the 
hyperfine coupling tensors A\ led to an angular- 
sensitive result. Inspired by this result, we 




FIG. 7: Entanglements for different angles under 
the hyperfine coupling tensors A\. 



propose that each electron interacts with ad- 
ditional local magnetic fields, Bi, rather than 
with the hyperfine fields. The Hamiltonian for 
this model is given by Eq. but with Av Ii 
replaced by Bi, the local magnetic field for the 
i th electron spin. We take the local fields to be, 
B x = (0,0, 4G), B 2 = (0,5G,0). 

In this case, we use the violation of the CHSH 
(Clauser, Home, Shimony, and Holt) inequal- 
ity [I?} as a witness of entanglement, a ver- 
sion of Bell's inequality 0. The CSHS in- 
equality is given by \E\ < 2X^ nax , where \E\ =| 
E(0,0) + E(0,t)+E(t,0)-E(t,t) |, E(t u t 2 ) = 
| • a)(f?2 • b) | <fit 2 ) is the two-time cor- 
relation function for a spin pair, and a and b 
are the unit direction vectors. The quantity, 
^max, is the maximum eigenvalue for the mea- 
surement operator, (<?i • a){B% • &), which for our 
specific operator is equal to 1. When \E\ ex- 
ceeds 2\^ lax = 2, the correlation between the 
two spins can no longer be explained classically, 
so the system is entangled. 

Fig. [5] shows the CHSH inequality as a func- 
tion of time for various orientations of the sys- 
tem in a magnetic field of 0.5G. Because there 
are now two perpendicular fields acting on the 
system, it becomes necessary to consider az- 



5 



imuthal orientation in addition to polar orien- 
tation. As seen in Fig. [8l as increases from 0° 
to 180°, the time for which the electron pair is 
entangled increases from roughly 60 ns to nearly 
90 ns, while for </> from 0° to 150° the variation 
of time of entanglement is restricted to an inter- 
val of less than 10 ns. It is interesting to note 
that this variation in time of entanglement oc- 
curs roughly on the same 100 ns time scale that 
the two electrons remain separated 

Changing the relative angles and strengths of 
the local magnetic fields has a dramatic impact 
on the angular sensitivity. A change in the field 
strength of the first electron from 4G to 5G is 
enough to dramatically increase both the az- 
imuthal and angular sensitivity of the entangle- 
ment. 

If indeed a protein such as cryptochrome is 
in part responsible for magnetoreccption, there 
must be some directional bias of the orientation 
of the protein, so that there will be a strong 
net signal. It is possible that this directional 
dependence could be provided by embedding 
within the membrane shelves of the photore- 
ceptor cells. This form of embedding leaves 
the protein free to rotate about one axis, but 
greatly restricts the rotation about its second 
axis [19|. For this reason it is necessary for the 
RP compass to be sensitive to rotation about 
one axis, while being virtually unaffected by ro- 
tation about the second. If the RP compass 
were to be sensitive to rotation in both 6 and (f>, 
the result of randomly oriented proteins about 
the axis would average out to create a back- 
ground signal that could potentially reduce the 
contrast of the RP compass. 

At the present time little is known about how 
cryptochrome is situated within the retina, in 
particular how it embeds onto or within the 
cell membrane [III]. There is no reason to as- 
sume that the z-axis of the RP model coincides 
with the fixed rotational axis of the embed- 
ded protein. As such, a configuration such as 
Si = (0, 0, 5G) and B 2 = (0, 5G, 0) might still 
produce a strong directional response under a 
coordinate transformation to the axis of protein 
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FIG. 8: Polar and azimuthal dependence of the 
CHSH inequality. The dashed blue line represents 
the points above which the system is entangled. 
The various orientations are given by the lines: 
red (0°), orange (30°), yellow (60°), green (90°), 
light blue (120°), dark blue (150°), black (180°). 
For the upper two figures, Bi = (0, 0, 5G) and 
B2 = (0, 5G, 0). The upper left figure depicts the 
azimuthal dependence for fixed 9 — 0°, while the 
upper right figure depicts polar dependence for fixed 
ip = 90°. For the lower two figures, Bi = (0,0, 5G) 
and £?2 = (0, 5G, 0). The bottom left and bottom 
right figures similarly depict the azimuthal and po- 
lar dependency, respectively. For all four figures 
k — 1/^s -1 . 



rotation. 

It should be pointed out that this model, un- 
like the previous model and its variants, is not 
symmetric about 90°, but it is symmetric about 
180°. While this might seem to contradict an 
inclination-only compass model, it is reasonable 
to assume that cryptochrome is cither bound to 
both sides of the cell membrane, or embedded 
within the membrane in both up and down ori- 
entations, so that the net signal cannot discern 
the polarity of the geomagnetic field. 

Conclusions and Future Work. — We have 
identified that the entanglement decay rate is 
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one of the key factors in the radical pair mecha- 
nism from the change of triplet yields (Fig. [4]). 
Also, we confirmed that the entanglement en- 
dures long enough for living systems to con- 
duct the entanglement-based reactions. How- 
ever, the dynamics of the entanglement is not 
sensitive to the change of angle between the z 
axis of the radical pair and the geomagnetic field 
vector in the hyperfine model. Therefore, if we 
still believe that entanglement plays a crucial 
role in the orientation of birds as demonstrated 
before, there must be an indirect mechanisms 
by which the entanglement can affect the birds' 
behavior. 

In the future, we will adjust the decay rates, 
for example, using different values for the decay 
rates of the singlet state and the triplet state 
to improve our model. We will also attempt 
to find the hidden bridge between the entangle- 
ment of the radical pair and the determination 
of orientation in a magnetic field. 
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